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摘 要 : 喀 喇 昆仑 山区 冰川 由 于 存在 正 物质 平衡 或 跃 动 .前进 现象 ,被 称 之 为 “ 喀 喇 昆仑 异常 ”, 不 过 
该 地 区 冰川 变化 差异 显著 ,尤其 是 大 型 表 磺 履 盖 冰川 ,呈现 与 其 他 类 型 冰川 明显 的 差异 性 响应 ,为 


理解 喀 咖 昆仑 冰川 异常 的 机 理 ,冰川 义 度 的 详细 变化 研究 十 分 人 必要 。 音 苏 盖 提 冰川 位 于 喀 喇 昆仑 


山 乔 艾 里 峰 北 坡 ,是 中 国 面 积 最 大 的 冰川 ,是 典型 的 大 型 表 磺 覆盖 冰川 。 通 过 应 用 TanDEM-X/Ter- 


raSAR-X(2014 年 2 月 ) 与 SRTM-X DEM(2000 年 2 月 ) 的 差分 干涉 测量 方法 计算 音 苏 盖 提 冰川 表 


面 高 程 变化 ,并 结合 冰川 表面 流速 对 冰川 表面 高 程 变化 和 路 动 进 行 分 析 和 讨论 。 结 果 表 明 :2000 一 
2014 年 音 苏 盖 提 冰 川 表面 高 程 平 均 下 降 了 1.68 €0.94 m, 即 冰川 整体 厚度 在 减 薄 , 年 变化 率 为 
-0.12+0.07 m .a 。 冰 川 表面 高 程 变化 分 布 不 均 , 其 中 南 分 支 (S) 冰 流 冰 川 整体 减 薄 较 为 显 


著 , 冰 川南 分 支 冰 流 运动 速度 较 快 ,前 进 / 跃 动 的 末端 占据 了 冰川 的 主干 , 阻 洁 原 主干 冰川 物质 的 向 
下 运 移 ( 跃 动 ) ,导致 原 主 干冰 舌 表面 高 程 上 升 ;冰川 厚度 减 薄 随 着 海拔 升 高 先 下 降 后 保持 稳定 , 同 
时 呈现 一 定 的 波动 性 ; 低 海 技 表 矿区 域 消融 大 于 裸 冰 区 ,可 能 存在 较 薄 表 矿 , 因 热传导 高 、 柳 盖 大 量 
冰 面 湖 塘 和 冰 崖 存在 ,加 速 了 冰川 消融 ;在 坡度 小 于 30。 的 区 域 ,冰川 厚度 减 薄 随 着 坡度 的 减 小 而 


加 剧 ; 坡 向 朝 南 冰川 厚度 略微 增加 (0.01 m) ,西南 坡 向 冰川 厚度 略微 减 薄 ( -0.03 m) ,其 他 坡 向 冰 
川 厚度 减 薄 明 显 。 近 14 a 来 , 表 奔 覆盖 的 音 苏 盖 提 冰川 表面 高 程 整体 下 降 表明 物质 处 于 亏损 状 


态 , 冰 川 跃 动 导 致 局 部 冰川 表面 高 程 的 增加 。 


X 键 词 : 音 苏 盖 提 冰 川 ; TanDEM-X/TerraSAR-X; 差分 干涉 测量 ; 表面 高 程 变化 ; 冰川 跃 动 


冰冻 圈 是 与 气候 变化 相互 影响 的 重要 圈 层 ， 
是 气候 变化 的 指示 器 "| ,也 是 重要 的 淡水 资源 宝 
Jg ,是 重要 的 服务 圈 层 ”; 。 气 候 变 暖 背景 下 ,全 球 
冰川 普遍 处 于 退缩 状态 ,而 喀 喇 昆仑 山大 部 分 冰川 
呈现 前 进 或 者 保持 稳定 , 即 “ 喀 喇 昆 仑 异常 。 
GARDELLE 等 利用 两 期 数字 高 程 模型 (DEM ) 进行 
研究 ,发 现 喀 喇 昆仑 山 冰 川 在 1999—2008 年 呈现 
正 物质 平衡 状态 ( +0.11+ 上 0.22 m w.e. a^!) 7, 
GARDNER 等 和 KAAB 等 分 别 利 用 ICESat GLAS 
数据 对 喀 喇 昆仑 山 冰川 表面 高 程 变化 进行 估算 ， 
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0.07m w. e. a^! ) ?' , fite ^ fii JH TanDEM-X 与 
SRTM dem 进行 差分 获得 结果 为 0.02 +0. 064 m - 
a^!" BOLCH 等 和 ZHOU 等 也 与 之 相当 结 
BE) 不 过 以 上 研究 结果 也 同时 表明 冰川 变 
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蒋 宗 立 等 :2000 一 2014 年 喀 喇 昆 


盖 冰 川 呈现 与 其 他 类 型 冰川 明显 的 差异 性 响应 ， 
表现 为 显著 的 物质 亏损 状态 。 表 研 尾 盖 冰 川 在 该 
区 域 分 布 广泛 ,与 裸 冰 相 比 , 表 奔 覆盖 层 由 于 反射 
率 .颗粒 大 小 .颜色 等 物理 性 质 的 差异 ,具有 独特 
的 热力 过 程 ,导致 其 下 履 冰 川 冰 的 消融 过 程 不 
同 “ ,使 得 表 研 覆 盖 冰川 与 其 他 类 型 冰川 对 同样 
的 气候 变化 条 件 显示 了 不 同 的 响应 特征 。 野 外 观 
测试 验 与 研究 表明 '"'"” , 表 奢 厚度 及 其 空间 分 布 
对 下 覆 冰 川 消融 过 程 的 影响 显著 ,从 而 影响 冰川 
物质 平衡 的 空间 分 布 特征 与 高 度 结构 特征 。 同 时 
大 型 表 磺 覆 盖 冰 川 冰 储 量 丰 富 ,对 区 域 河川 径流 
和 水 资源 有 重要 影响 ,在 气候 变 暧 的 影响 下 ,冰川 
储量 的 变化 影响 着 下 游 水 资源 和 灾害 风险 。 此 外 
喀 喇 昆仑 山 的 冰川 跃 动 也 存在 不 同 的 机 理 o 
此 选择 大 型 表 奢 覆盖 型 冰川 进行 详细 人 研究 十 分 必 
要 。 音 苏 盖 提 冰川 位 于 喀 喇 昆仑 山中 部 ,是 中 国 
境内 覆盖 面积 最 大 的 冰川 ,也 是 典型 的 大 型 表 硬 
覆盖 冰川 ,该 冰川 分 支 发 生 过 多 次 路 动 , 据 COP- 
LAND 等 的 研究 表明 ,该 冰川 至 少 有 8 处 发 生 过 
跃 动 或 者 疑似 跃 动 ,其 中 5 处 确认 发 生 过 跃 动 。 
蒋 宗 立 等 和 刘 星 月 等 针对 音 苏 盖 提 冰川 运动 速度 
进行 了 详细 研究 “~ ,但 目前 还 缺乏 与 音 苏 盖 提 
冰川 高 程 变化 的 相 结合 研究 ,集成 研究 结果 将 有 
助 于 提高 对 大 型 表 磺 覆盖 冰川 的 冰川 变化 差异 响 
应 和 冰川 路 动机 制 的 认识 。 

由 于 音 苏 盖 提 冰川 深 居 克 勒 青河 中 上 游 ,地 


面 补 给 困难 、 人 力 难 以 企及 ,而 遥感 技术 是 目前 切 
实 可 行 的 研究 手段 。 近 些 年 随 着 光学 遥感 技术 的 
不 断 发 展 ,在 高 分 辨 率 DEM 数据 提取 方面 取得 了 
较 大 进展 ,但 是 受 云 遮 等 天 气 条 件 所 限制 , 且 冰 和 雪 
对 于 可 见 光 的 反射 率 过 高 ,难以 在 冰雪 覆盖 区 域 
获取 特征 点 ,获取 冰雪 区 数字 高 程 模型 误差 很 大 。 
合成 孔径 雷达 (Synthetic Aperture Radar,SAR ) 主 
动 微波 传感器 不 受 云 雪 影 响 ,全 天 候 全 天 时 工作 
可 弥补 光学 影像 的 不 足 , 国 内 外 学 者 利用 合成 孔 
径 雷 达 干 涉 测量 (InSAR ) 方 法 获取 了 大 范围 高 精 
度 的 数字 高 程 模 型 数据 ,并 且 广 泛 应 用 于 冰川 变 
AL ERS 0777. AN IC All J] TanDEM-X/TerraSAR- 
X( TDX/TSX ) 双星 系统 的 SAR 数据 和 SRTM-X 
DEM 数据 ,采用 D-InSAR 方法 获取 了 2000—2014 
年 音 苏 盖 提 冰川 表面 高 程 变化 信息 ,并 对 其 冰川 
表面 高 程 变 化 特征 进行 分 析 讨 论 。 
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仑 山 音 苏 盖 提 冰川 表面 高 程 变化 


1 研究 区 概况 


喀 喇 昆仑 山 是 全 球 中 低 纬度 地 区 山地 冰川 分 布 
数量 最 多 ,冰川 规模 最 大 的 区 域 之 一 。 该 区 域 位 于 
亚 欧 大 陆 中 心 , 属 于 典型 大 陆 型 气候 ,气温 年 较 差 
大 ,降水 稀少 。 但 由 于 该 区 域 位 处 西风 带 ,又 受到 减 
弱 印 度 洋 西 南 季风 的 影响 ,有 冰川 分 布 的 高 海拔 地 
区 年 降水 量 可 达 2 000 mm 以 上 。 音 苏 盖 提 冰川 (中 
心 位 置 坐标 为 36.083"N ,76. 1028°E) 位 于 喀 喇 昆仑 
山中 部 乔 苞 里 峰 西 侧 北 坡 .克勤 青河 中 段 的 南岸 ,总 
面积 为 493 km? ,是 我 国 面积 最 大 的 冰川 55 (图 1)。 
冰川 主体 呈 自 西向 东 的 流向 ,由 4 条 主要 冰 流 (SW、 
NW N.S) 和 超过 10 条 小 分 支 冰 流 汇流 而 成 ,冰川 
主流 线 最 大 长 度 达 42 km ,末端 海拔 约 4 000 m, 西 
南 分 支 冰 流 (SW) 为 冰川 主 源头 (图 1)。 
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注 :N 为 北 冰 流 ,NW 为 西北 冰 流 ,SW 为 西南 冰 流 ,S 为 南 冰 流 ,已 报 


道 的 跃 动 冰川 序号 与 Copland 等 研究 保持 一 致 
图 1 音 苏 盖 提 冰川 位 置 范围 


Fig.1 Location and extent of Yengisogat Glacier 
2 数据 来 源 与 研究 方法 


2.1 数据 来 源 

本 文 主要 采用 航天 飞机 地 形 测绘 任务 (Shuttle 
Radar Topographic Mission ,SRTM ,2000 年 2 月 )X- 波 
Et DEM 与 TanDEM-X/TerraSAR-X ( TDX/TSX ) 数据 
(2014 年 2 月 ) 对 冰川 表面 高 程 变化 进行 分 析 ( 表 
1)。SRTM 是 美国 宇航 局 (NASA ) 2000 年 发 射 的 
“奋进 "号 航天 飞机 上 搭载 的 对 地 观测 雷达 系统 ,在 
轨 高 度 233 km ,使 用 C 波段 和 X 波段 分 别 采用 225 
km 和 50 km 的 幅 宽 进行 地 面 扫描 ,获取 了 56°S 到 
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tió wni 
R1 数据 源 7000 r 
Tab.1 Information of data sources used J-LO hia? 
6500 上 
数据 获取 时 间 工作 波段 “分 辩 率 / m © 
TDX/TSX 2014 年 2 月 8 日 X 波段 8 x10 g 6000 f 
2014 4£2 H 19 H D 
| & 5500 上 
SRTM DEM 2000 年 2 H C 波段 30 x30 E 
XE & 5000 + 
o MN B3 Y hy at > 3 Y 4500 f 
60°N 之 间 的 DEM 数据 ,空间 分 辩 率 为 1 弧 秒 , 约 为 
30 m"* , TDX/TSX 系统 是 德国 宇航 局 (DLR ) 2007 4 000 | 
4000 4500 5000 5500 6000 6500 7000 


年 6 月 和 2010 年 6 月 分 别 发 射 的 两 颗 SAR 卫星 组 
成 的 单轨 双 天 线 双星 分 布 式 系统 ,在 轨 高 度 514. 8 
km ,重复 周期 11 d, 聚 束 模式 下 分 辨 率 可 达 E m”, 
本 研究 使 用 该 系统 双 基 站 strip map 模式 获取 的 精 
确 配 准 的 单 视 复 (Single Look Complex , SLC ) 数据 
对 。 
2.2 ”研究 方法 

使 用 1 弧 秒 SRTM dem ,像素 大 小 为 约 30 m ,对 
TanDEM-X/TerraSAR-X 干涉 图 进行 12 x15 多 视 处 
理 ,与 SRTM dem 像素 大 小 一 致 。PmSAR 生成 的 干 
涉 相位 中 包括 地 形 相 位 、 基 线 误 差 导致 的 相位 JE AE 
相位 ,大气 延 迟 引 起 的 相位 及 噪声 相位 , 因 双 基站 模 


式 的 大 气 延迟 .噪声 相似 ,干涉 计算 之 后 可 以 忽略 ， 
也 不 存在 形变 相位 。 为 得 到 地 表 高 程 变化 ,通过 差 
分 处 理 , 去 除 SRTM DEM 模拟 的 地 形 相 位 ,并 去 除 
基线 误差 导致 的 线性 趋势 相位 ,从 而 得 到 冰川 表面 
高 程 变化 的 残 差 相位 ,把 残 差 相位 转化 为 高 程 变化 
值 , 即 可 获得 冰川 表面 的 变化 |。 

本 文采 用 的 参考 DEM 是 2000 年 获取 的 SRTM- 
X dem 数据 ,相对 于 SRTM C dem, SRTM-X dem Il 
声 比较 大 ; 因 SRTM-X 或 C DEM 在 喀 喇 昆仑 山坡 度 
很 大 的 区 域 存在 大 量 数据 空洞 ,这 些 空 值 区 域 分 布 
在 海拔 较 高 区 域 , 即 冰川 积累 区 ,予以 剔除 。 但 在 数 
据 预 处 理 时 使 用 填补 了 空 值 区 域 的 完整 覆盖 数字 高 
程 模 型 。 

由 于 两 对 TSXZTDX 影像 获取 时 间 为 2 月 份 ,而 
且 仅 仅 相 差 11 d, 本 文 忽略 此 时 间 段 的 冰川 高 程 变 
化 。 不 过 由 于 轨道 偏 移 等 因素 的 影响 ,两 对 影像 的 
计算 结果 在 重合 区 域 有 一 定 误差 ,为 了 消除 该 部 分 
误差 ,本文 对 两 幅 影 像 结 果 的 重 县 区域 误 差 进 行 了 
统计 分 析 , 结 果 发 现 其 间 存 在 显著 的 线性 关系 (图 
2)。 因 此 根据 拟 合 得 到 的 线性 关系 式 以 2014 年 2 
月 8 日 的 数据 结果 为 基准 对 2014 年 2 H 19 日 的 数 
据 结 果 进 行 高 程 调整 , 使 得 两 幅 结 果 图 可 以 无 缝 
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图 2 EDS BRU 


Fig.2 Fitted trend line of elevation error 


拼接 。 此 外 , 选 在 非 冰 川 相 对 平坦 区 域 进行 计算 结 
果 的 误差 计算 "” ,最 终 误差 为 : £0.07 m a! 


3 结果 与 分 析 


2000—2014 年 音 苏 盖 提 冰川 表面 高 程 平均 下 
降 了 1. 68 +0.94 m, 即 冰川 整体 厚度 在 减 薄 ,年 变 
化 率 为 -0.12 £0.07 m .a ”。 冰 川 表面 高 程 变 化 
分 布 不 均 ,其 中 北 冰 流 CN) 上游 (图 3a 和 图 4a) \ 西 
北 冰 流 (NW) 上 游 ( 图 3b 和 图 4b) 北 冰 流 中 部 (图 
3c 和 图 4a) 西北 冰 流 中 部 (图 3d 和 图 4b) 和 干流 
中 部 (西南 冰 流 ,图 3e 和 图 4c) 均 有 明显 的 冰川 增 
厚 现象 , 其 中 干流 中 部 增 厚 最 为 显著 ,最 大 增 厚 
45.77 m。 南 冰 流 冰川 (图 4d) 减 薄 最 为 显著 ,其 中 
末端 变化 尤为 强烈 。 
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图 3 2000—2014 年 冰川 表面 高 程 变化 
Fig.3 Results of elevation changes between 2000 and 2014 
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图 4 音 苏 盖 提 冰川 主要 分 支 表面 高 程 变化 剖面 图 


Fig.4 Surface elevation change profiles along the main flows of Yengisogat Glacier 


4 讨论 


4.1 冰川 高 程 变化 影响 因素 


海拔 4 700 m 以 上 到 雪线 的 冰 体 发 生 印 载 转移 至 海 
Hk 4 500 ~4 700 m 之 间 。 HEH 4 500 m 以 下 冰川 厚 
度 减 薄 随 着 海拔 升 高 较为 显著 , 受 冰 川 运 动 的 影响 ， 


降水 与 气温 是 影响 冰川 发 育 的 主要 气候 因子 ， 
温度 影响 冰川 消融 ,降水 影响 冰川 积累 ,两 者 变化 共 
同 决 定 冰 川 的 性 质 \ 发 育 与 变化 。 周 围 气象 站 资料 
表明 , 近 些 年 该 地 区 气温 和 降水 都 有 升 高 的 趋势 , 气 
温 升 高 导致 冰川 消融 ,降水 增加 又 在 一 定 程度 上 弥 
补 了 冰川 的 消融 。 由 于 喀 喇 昆仑 山 气象 观测 资料 比 
较 稀 少 , 并 且 局 部 气候 较为 复杂 ,无 法 对 冰川 高 程 变 
化 与 气候 变化 进行 定量 分 析 。 不 过 一 般 情况 下 , 气 
温 变 化 对 冰川 物质 平衡 起 主导 作用 ,一 方面 气温 每 
上 升 1 的 物质 亏损 需要 降水 增加 25% 才能 弥补 ， 
另 一 方面 气温 升 高 导致 液态 降水 比例 增多 以 及 加 快 
冰川 表面 的 融雪 过 程 ,降低 冰川 表面 反照 率 , 促 进 冰 
川 消融 。 本 研究 中 冰川 表面 高 程 降 低 , 主 要 是 受气 
温 升 高 所 影响 。 

地 形 影 响 水 热 分 布 ,从 而 影响 冰川 的 积累 和 消 
融 。 研 究 结果 表明 (图 5) , 音 苏 盖 提 冰 川 雪 线 海拔 
高 度 在 5 310 ~5 500 m zz [i] ,雪线 海拔 以 上 表面 高 
程 基本 保持 不 变 。 雪 线 以 下 海拔 4 500 ~ 4 700 m 之 
间 冰 川 表面 明显 增 厚 ,此 部 分 增 厚 是 由 跃 动 造成 的 ， 


高 程 变化 的 海拔 依赖 关系 也 受到 一 定 的 波动 。 低 海 
拔 表 磺 区 域 的 冰川 减 薄 比 裸 冰 区 显著 ,可 能 是 此 区 
域 表 磺 相对 较 注 , 薄 表 研 具 有 较 小 的 热 容 ,内 部 具有 
很 高 的 热传导 率 ,从 而 促进 冰川 消融 ,同时 薄 表 研 区 
域 冰 面 湖 、 冰 崖 等 也 较为 发 育 ,伴随 着 热 喀斯特 作用 
对 冰 面 具有 加 剧 消融 作用 。 

坡度 小 于 30° 时 ,冰川 厚度 减 薄 随 着 坡度 的 减 
小 而 加 剧 ,较为 平坦 的 区 域 一 般 海拔 较 低 , 且 坡 度 较 
组 区域 冰川 运动 速度 较 慢 ,受到 高 海拔 地 区 的 补给 
较 少 (图 6)。 坡 度 为 30 ° ~50 "时 ,冰川 厚度 变化 较 
小 ,此 区 域 一 般 在 积累 区 ,积累 区 变化 较 小 。 当 坡度 
大 于 50 "时 ,冰川 厚度 随 着 坡度 的 增加 而 加 速 减 薄 ， 
该 区 域 坡度 陡峭 , 易 发 生 雪 骨 或 者 冰 骨 ,从 而 导致 表 
面 高 程 减 薄 , 不 过 该 区 域 冰川 面积 非常 小 , 较 陡峭 的 
区 域 本 身 就 不 利于 雪 冰 的 积累 。 

由 图 7 可 知 , 南 坡 冰 川 厚度 略微 增加 (0. 01 
m) ,西南 坡 向 冰川 厚度 略微 减 薄 ( -0.03 m) ,其 他 
坡 向 冰川 厚度 减 薄 明显 ,尤其 是 北 坡 和 西北 坡 。 南 
坡 和 西南 坡 主要 对 应 北 冰 流 和 西北 冰 流 , 北 坡 主要 
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图 7 音 苏 盖 提 冰川 不 同 朝向 的 平均 高 程 变化 


Fig.7 Average elevation change for different aspect of 


Yengisogat Clacier 


对 应 朝 南 冰 流 和 朝 西 南 冰 流 , 从 图 3 空间 分 布 格局 
上 也 可 以 看 出 , 朝 北 冰 流 和 西北 冰 流 积累 区 明显 增 
JE ,冰川 表面 高 程 整体 变化 较 小 ,这 是 朝 南 冰 流 接收 
更 多 太阳 辐射 的 缘故 。 

北 冰 流 上 游 和 西北 冰 流 上 游 均 位 于 或 者 接近 积 
累 区 , 陡 壁 下 山谷 内 的 冰川 增 厚 尤为 明显 ,该 部 分 冰 
川 物质 很 可 能 来 源 于 陡 壁 处 的 雪 月 和 各 小 分 支 冰川 
的 汇流 。 西 北 冰 流 中 部 、 北 冰 流 中 部 和 干流 中 部 的 
上 端 冰川 有 明显 的 减 薄 ,因此 可 认为 其 物质 分 别 大 
幅度 转移 至 西北 冰 流 中 部 、 北 冰 流 中 部 和 干流 中 部 ， 
其 物质 转移 可 能 是 受 跃 动 所 驱动 。 也 可 能 是 南 分 支 
冰 流 几乎 阻 断 了 主干 冰 流 ,冰川 物质 在 此 堆积 , 厚 表 
磺 覆 盖 使 得 消融 减弱 ,冰川 厚度 增加 ,该 区 域 冰 川 厚 
度 的 持续 增加 ,有 可 能 导致 主干 向 前 推进 ,重新 成 为 
主 冰 流 。 南 冰 流 冰川 (图 6b ) 减 薄 最 为 显著 ,其 中 末 
端 变 化 尤为 强烈 。 南 冰 流 ( 南 斯 嘎 姆 里 冰川 ) 在 
1978—1990 年 间 发 生 过 跃 动 ,并 取代 西南 冰 流 的 主 
导 地 位 。 虽 然 南 冰 流 冰 川 不 在 跃 动 期 ,但 冰川 流速 
较 大 ,COPLAND 等 .1 研究 发 现 2006—2007 年 南 冰 
流 消融 区 上 游 冰 川 运动 速度 大 于 200 m a”, 且 末 
流速 度 逐 渐 减 小 ,将 宗 立 等 和 刘 星 月 等 的 研究 结 
也 与 之 类 似 “-”。 较 快 的 冰川 流速 导致 冰川 物质 
不 断 由 上 游 转移 到 下 游 , 低 海拔 冰川 更 容易 受到 和 气 
温 升 高 的 影响 ,因此 南 冰 流 冰川 整体 减 薄 较为 显著 。 
与 此 同时 , 南 分 支 以 下 区 域 至 冰川 末端 冰 流 速度 却 
非常 微弱 ,导致 末端 冰川 得 不 到 更 多 的 物质 补给 。 
因此 ,在 气温 变 暖 的 影响 下 ,末端 及 附近 区 域 冰 川 减 
薄 最 为 显著 。 结 合 前 期 获得 的 流速 分 布 特征 ,我们 
推断 :有 可 能 是 南 分 支 冰 流 几乎 阻 断 了 主干 冰 流 ,上 
游 传输 的 冰川 物质 在 此 堆积 ,加 上 厚 表 矿 履 盖 使 得 
消融 减弱 ,冰川 厚度 增加 ,该 区 段 冰川 厚度 的 持续 增 
加 ,有 可 能 导致 主干 向 前 推进 ,重新 成 为 主 冰 流 。 
4.2 冰川 跃 动 

冰川 跃 动 通常 是 指 冰 川 周期 性 地 在 较 短 时 间 内 
运动 速度 突然 加 快 并 引起 冰川 物质 快速 转移 和 冰 体 
重新 分 布 的 现象 。 由 图 3 可 以 看 出 ,c、d、e 处 冰 体 
明显 增 厚 ,同时 上 游 出 现 了 明显 减 薄 ,符合 冰川 路 动 
村 征 ,其 中 e 处 最 为 明显 。 将 宗 立 等 研究 发 现 西南 
冰 流 在 2007—2009 年 冰川 运动 速度 很 快 (100 m - 
a^) ,QUINCEY 等 认 结 果 表 明 该 冰 流 2012 年 6~8 
月 运动 速度 约 900 m * a^! ,随后 冰川 运动 速度 快速 
下 降 , 即 2012 年 6 ~8 月 可 能 处 于 跃 动 期 最 高 峰 。 
COPLAND 等 认为 该 冰川 1973 年 以 前 发 生 过 跃 动 ， 
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本 研究 无 法 判断 出 该 冰川 在 此 发 生 跃 动 的 具体 时 
间 , 不 过 根据 现 有 资料 ,2012 年 发 生 跃 动 的 可 能 性 
大 一 些 。 图 3 的 c 人 处 和 d 处 跃 动 特征 不 显著 ,c 处 
此 前 没有 报道 过 跃 动 现象 ,不 过 其 所 在 北 音 苏 盖 
提 冰 川 有 多 次 跃 动 或 疑似 跃 动 现象 ,其 中 c 处 西 
南 侧 文 (图 1 中 26 号 点 )1977 一 1990 年 曾 发 生 过 
ERZ) sd 处 此 前 也 没有 跃 动 现象 的 报道 ,其 西南 侧 
文 1990—2001 F Ñ RETR KWR T 
COPLAND 等 研究 结果 ,认为 音 苏 盖 提 冰川 至 少 有 
11 处 跃 动 或 者 疑似 跃 动 迹象 ,COPLAND 等 的 结 
确认 了 5 Ra ^ ,本 研究 确认 西南 冰 流 (图 3e 
处 ) 为 跃 动 。 

目前 ,冰川 跃 动 的 原因 主要 有 两 种 机 制 ,一 是 热 
力学 机 制 (多 温 型 冰川 ,比如 斯 瓦尔 巴特 群岛 类 型 
的 冰川 ) ,三 是 水 文学 机 制 ( 温 冰 川 ,比如 阿拉 斯 加 
SFY AVI] 7; QUINCEY 等 "7” 认为 喀 喇 昆仑 山 
的 冰川 不 受 任何 单一 机 制 所 控制 ,并 且 与 水 文 和 地 
貌 条 件 有 关 。 对 于 高 差 极 大 的 冰川 来 说 ,其 补给 区 
位 于 极 高 海拔 ,在 消融 区 可 下 降 到 很 低 的 海拔 高 度 ， 
其 冰川 温度 的 分 布 类 型 变化 很 大 ,往往 在 冰川 上 游 
为 冷 型 ,而 在 冰川 下 游 为 温 型 , 整 条 冰川 应 该 属于 多 
温 型 ,比如 喀 喇 昆仑 山 的 巴 托 拉 冰 川 *, 音 苏 盖 提 
冰川 也 可 能 属于 多 温 型 冰川 。 音 苏 盖 提 冰川 是 典型 
的 树枝 状 冰川 ,其 各 个 分 支 的 性 质 可 能 存在 差异 AR 
据 QUINCEY 等 "的 研究 结果 ,西南 冰 流 跃 动 峰值 
速度 较 大 ,2012 年 夏季 运动 速度 达到 峰值 ,然后 运 
动 速度 趋 于 正常 , 跃 动 阶 时 间 较 短 ( 几 个 月 ), 不 太 
符合 热力 机 制 引起 的 跃 动 。 值 得 一 提 的 是 , 音 苏 盖 
提 冰 川 被 称 为 裂隙 冰川 , 冰 面 裂隙 广泛 分 布 ,如 果 液 
态 降水 或 者 冰川 融 水 从 裂隙 进入 冰 体 ,有 可 能 改变 
其 水 文 条 件 ,水 的 润滑 促进 了 冰川 跃 动 。 
4.3 与 其 他 研究 的 比较 

目前 针对 喀 喇 昆仑 山 冰 川 表面 高 程 变化 或 者 物 
质 平衡 研究 中 ,普遍 认为 喀 喇 昆仑 山 冰 川 整体 上 呈 
现 接近 于 平衡 的 状态 , 不 过 空间 差异 性 很 大 。 
BRUN 等 利用 2000 一 2016 年 ASTER DEM 得 到 音 苏 
盖 提 冰川 高 程 变化 趋势 为 -0.11 ma ,冰川 减 薄 
速率 高 于 喀 喇 昆仑 山 2000—2016 年 平均 值 -0. 03 
m.a hl ,与 本 研究 结果 ( -0.12 m + a7!) — 3K, 
ZHOU 等 (1 利用 KH9 和 SRTM 的 研究 结果 表明 音 
苏 盖 提 冰川 所 在 的 喀 喇 昆仑 山中 东部 1979—2000 
年 -0.19+ 上 0.08 m* a ,高 于 整个 喀 喇 昆仑 山 的 
1979—2000 年 平均 值 ( -0.09 £0.03 m.a-)。 


仑 山 音 苏 盖 提 冰川 表面 高 程 变化 


此 , 音 苏 盖 提 所 在 区 域 冰 川 亏 损 比 咯 喇 昆 仑 山 其 他 
区 域 高 ,临近 的 一 些 冰 川 其 至 呈正 平衡 状态 ,比如 
Baltoro 冰川 1999—2010 年 冰川 物质 平衡 为 +0. 10 
+0.13 mw.e. a^! ,该 冰川 在 此 期 间 没 有 跃 动 的 
迹象 ,除了 地 形 差异 外 , 跃 动 后 冰川 物质 从 上 游 转移 
到 下 游 , 低 海拔 冰 体 受气 温 影 响 更 敏感 ,物质 亏损 相 
对 较 大 , 音 苏 盖 提 冰川 跃 动 比较 频繁 , 跃 动 后 冰川 表 
面 高 程 减 薄 速度 加 快 。 


5 结论 


鉴于 目前 缺乏 冰川 尺度 的 大 型 表 研 覆盖 型 冰川 
高 程 变化 的 详细 研究 ,本 文 基于 TanDEM-X/Terra- 
SAR-X (2014 年 2 月 ) 与 SRTM-X DEM (2000 4E 2 
月 ) 的 差分 干涉 测量 方法 进行 了 音 苏 盖 提 冰 川 表面 
高 程 变化 的 研究 ,结果 表明 : 

(1) 2000—2014 年 音 苏 盖 提 冰川 表面 高 程 平 
均 下 降 了 1.68 +0.94 m, 即 冰川 整体 厚度 在 减 薄 ， 
年 变化 率 为 -0.12 +0.07 ma”。 冰 川 表面 高 程 
变化 分 布 不 均 , 其 中 南 冰 流 冰川 整体 减 薄 较为 显著 ， 
南 冰 流 冰 川 运动 速度 较 快 ,高 海拔 的 冰川 物质 转移 
到 低 海 拔 区 域 ,加 快 了 消融 的 进程 。 

(2) 气温 升 高 是 冰川 减 蒲 的 主因 ,地 形 影 响 水 
热 条 件 ,进而 影响 冰川 厚度 变化 。 冰 川 厚度 减 薄 随 
着 海拔 升 高 先 下 降 后 保持 稳定 ,但 也 呈现 一 定 的 波 
动 性 。 低 海拔 表 磺 区 域 消融 大 于 裸 冰 区 , 注 表 磺 热 
容 小 .热传导 高 AEE ait eek TEL TH BA Ke, PEE T 
冰川 消融 。 坡 度 小 于 30° 时 ,冰川 厚度 减 薄 随 着 坡 
度 的 减 小 而 加 剧 。 南 坡 和 西南 坡 向 冰川 厚度 几乎 不 
变 , 其 他 坡 向 冰川 厚度 减 薄 明 显 。 

(3) 本 研究 补充 了 COPLAND 等 音 苏 盖 提 冰川 
跃 动 事件 ,认为 2000—2014 年 音 苏 盖 提 冰川 有 3 处 
跃 动 或 疑似 跃 动 的 现象 。 

即使 喀 喇 昆仑 山 冰 川 整体 微弱 物质 增加 ,对 于 
表 磺 履 盖 的 音 苏 盖 提 冰川 而 言 , 近 14 a 来 处 于 物质 
亏损 状态 。 


致谢 :感谢 DLR AO 科学 用 户 (syl_cas_XTI_ 


LAND6642 ) 提供 TanDEM-X/TerraSAR-X 干涉 数据 
与 SRTM-X dem 数据 。 


参考 文献 ( References ) 


[1] OERLEMANS J, REICHERT B K. Relating glacier mass balance to 


meteorological data by using a seasonal sensitivity characteristic 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[14] 


td £4 


[J]. Journal of Glaciology ,2000 ,46(152) :1 -6. 

件 建 新 , 李 忠 勤 , 张 慧 ,等 . 中 国 西部 大 陆 性 冰川 与 海洋 性 冰川 
物质 平衡 变化 及 其 对 气候 响应 一 一 以 乌 源 1 号 冰川 和 帕 隆 94 
号 冰川 为 例 [ 可 .干旱 区 地 理 ,2019,42 (1):20 - 28. [MU 
Jianxin ,LI Zhongqin , ZHANG Hui,et al. Mass balance variation of 


continental glacier and temperature glacier and their response to 
climate change in western China; Taking Urumuqi Glacier NO. 1 
and Parlung NO. 94 Glacier as examples[ J]. Arid Land Geogra- 
phy 2019 ,42(1) :20 -28. ] 

Xi idi , AJI, PRIE BBE. 冰冻 圈 水 资源 服务 功能 区 划 研 究 综 述 
[J]. FA dh BE, 2018 ,41 (4) :751 -760. [LIU JIE, HUANG 


Jinchuan , LIN Haoxi. A literature review on regionalization of water 
resource service function of cryosphere in China[ J]. Arid Land 
Geography ,2018 ,41 (4) :751 -760. ] 

HEWITT KENNETH. The Karakoram anomaly? Glacier expansion 
and the ‘ elevation effect’ , Karakoram Himalaya[ J]. Mountain Re- 
search and Development ,2005 ,25 (4) :332 - 340. 

GARDELLE J, BERTHIER E, ARNAUD Y. Slight mass gain of 
Karakoram glaciers in the early twenty-first century[ J ] . Nature ge- 
oscience ,2012 ,5(5) :322 - 325. 

COPLAND L,SYLVESTRE T, BISHOP M P,et al. Expanded and 
recently increased glacier surging in the Karakoram [J]. Arctic, 
Antarctic ,and Alpine Research ,2011 ,43 (4) :503 —516. 
GARDNER A S, MOHOLDT G,COGLEY J G,et al. A reconciled 
estimate of glacier contributions to sea level rise; 2003 to 2009 
[J]. Science ,2013 ,340 (6137 ) :852 - 857. 

KAAB A,TREICHLER D,NUTH C ,et al. Brief communication ; 
Contending estimates of 2003—2008 glacier mass balance over the 
Pamir-Karakoram-Himalaya[ J]. The Cryosphere ,2015 ,9 (2) :557 
- 564. 

BRUN F, BERTHIER E, PATRICKW, et al. A spatially resolved 
estimate of High Mountain Asia glacier mass balances from 2000 to 
2016[ J]. Nature Geoscience ,2017 ,10(9) :668 - 673. 

LIN H,LI G,CUO L,et al. A decreasing glacier mass balance gra- 
dient from the edge of the upper Tarim Basin to the Karakoram 
during 2000—2014[ J]. Scientific Reports ,2017 ,7 :6712. 

BOLCH T,PIECZONKA T, MUKHERJEE K ,et al. Brief communi- 
cation; Glaciers in the Hunza catchment ( Karakoram) have been 
nearly in balance since the 1970s [ J]. The Cryosphere, 2017 , 11 
(1) :531 -539. 

ZHOU Y S, LI Z W,LI J, et al. Geodetic glacier mass balance 
(1975—1999) in the central Pamir using the SRTM DEM and 
KH-9 imagery[ J]. Journal of Glaciology ,2019 ,65 :309 - 320. 
GARDELLE J, BERTHIER E, ARNAUD Y, et al. Region-wide 
glacier mass balances over the Pamir-Karakoram-Himalaya during 
1999—2011[ J]. The Cryosphere ,2013 ,7 (4) :1263 — 1286. 
ÓSTREM G. Ice melting under a thin layer of moraine , and the ex- 
istence of ice cores in moraine ridges[ J ]. Geografiska Annaler, 
1959 ,41 (4) :228 -230. 

BENN D I, BOLCH T, HANDS K, et al. Response of debris-cov- 


ered glaciers in the Mount Everest region to recent warming, and 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[27] 


[28] 


ChinaXiv& ERAT! 


implications for outburst flood hazards| J]. Earth-Science Reviews, 
2012,114(1 -2) :156 174. 

ZHANG Y,FUJITA K,LIU S Y ,et al. Distribution of debris thick- 
ness and its effect on ice melt at Hailuogou glacier, southeastern 
Tibetan Plateau, using in situ surveys and ASTER imagery [ J]. 
Journal of Glaciology ,2011 ,57(206) :1147 —1157. 

QUINCEY D J, GLASSER N F,COOK S J,et al. Heterogeneity in 
Karakoram glacier surges [ J ]. Journal of Geophysical Research; 
Earth Surface ,2015 ,120(7) :1288 - 1300. 

JIANG Z L,LIU S Y,PETERS J,et al. Analyzing Yengisogat Glac- 
ier surface velocities with ALOS PALSAR data feature tracking, 
Karakoram, China [ J ]. Environmetal Earth Sciences, 2012 , 67 
(4) 11033 - 1043. 


刘 星 月 , 蒋 宗 立 , 刘 时 银 y SE. 音 苏 盖 提 冰川 表面 流速 特征 分 析 
[J]. 地 球 学 报 ,2018 ,39(1) :103 — 110. [ LIU Xingyue, JIANG 


Zongli, LIU Shiyin, et al. An analysis of surface velocities on 
Yengisogat glacier, Karakoram China[ J]. Acta Geoscientica Sini- 
ca,2018 ,39(1) :103 -110] 

ZEBKER H A,GOLDSTEIN R M. Topographic mapping from in- 
terferometric synthetic aperture radar observations [ J]. Journal of 
Geophysical Research ,1986 ,91 ( B5) :4993 — 4999. 

KWOK R,FAHNESTOCK M A. Ice sheet motion and topography 
from radar interferometry [ J ]. IEEE Transactions on Geoscience 
and Remote Sensing , 1996 ,34(1) :189 —200. 

RANKL M, BRAUN M. Glacier elevation and mass changes over 
the central Karakoram region estimated from TanDEM-X and 
SRTM/X-SAR digital elevation models[ J]. Annals of Glaciology , 
2016 ,51(71) :273 -281. 

RIZZOLI P, MARTONE M, GONZALEZ C, et al. Generation and 
performance assessment of the global TanDEM-X digital elevation 
model[ J]. ISPRS Journal of Photogrammetry and Remote Sensing, 
2017,132:119 - 139. 

LAMBRECHT A,MAYER C, WENDT A, et al. Elevation change 
of Fedchenko Glacier, Pamir Mountains , from GNSS field measure- 
ments and TanDEM-X elevation models , with a focus on the upper 
glacier[ J]. Journal of Glaciology ,2018 ,64 (246) 637 - 648. 

刘 时 银 , 姚 晓 军 , 郭 万 钦 , 等 . 基于 第 二 次 冰川 编目 的 中 国 冰川 
BUR CJ]. 地 理学 报 ,2015 ,70(1) :3 -16. [ LIU Shiyin, YAO Xi- 


aojun, GUO Wanqin, et al. The contemporary glaciers in China 


based on the Second Chinese Glacier Inventory [ J ]. Acta Geo- 
graphica Sinica ,2015 ,70(1) :3 -16. ] 

MARSCHALK U, ROTH A, EINEDER M, et al. Comparison of 
DEMs derived from SRTM/X and C-Band [ C ] // Conference ; The 
Shuttle Radar Topography Mission-Data Validation and Applica- 
tions , Reston , Virginia , USA ,2005 , Cite this publication. 

RIZZOLI P, BRAUTIGAM B, KRAUS T, et al. Relative height er- 
ror analysis of TanDEM-X elevation data [J]. ISPRS Journal of 
Photogrammetry and Remote Sensing,2012 ,73 :30 - 38. 

NECKEI N, BRAUN A, KROPACek J,et al. Recent mass balance 
of the Purogangri Ice Cap, central Tibetan Plateau ,by means of dif- 


ferential X-band SAR interferometry [ J |. The Cryosphere ,2013 ,7 


蒋 宗 立 等 :2000 一 2014 年 喀 喇 昆仑 山 音 苏 盖 提 冰川 表面 高 程 变 化 


(5) :1623 - 1633. Svalbard and other regions [ J]. Journal of Geophysical Research, 
[29] COPLAND L, POPE S, BISHOP M P, et al. Glacier velocities 2003 , 108 ( B5) ;2237 ,doi:10. 1029/2002JB001906. 

across the central Karakoram[ J]. Annals of Glaciology , 2009 , 50 [31] 谢 自 楚 , 刘 朝 海 . 冰川 学 导论 LM]. 上 海 : 上 海 科学 普及 出 版 

(52) :41 -49. 社 ,2010:490. [ XIE Zichu,LIU Chaohai. Introduction of Glaciolo- 
[30] MURRAY T, STROZZI T, LUCKMAN A, et al. Is there a single gy[ M]. Shanghai; Shanghai Popular Science Publishing House, 

surge mechanism? Contrasts in dynamics between glacier surges in 2010 :490. ] 


Surface elevation changes of Yengisogat Glacier between 2000 and 2014 


JIANG Zong-li', WANG Lei', ZHANG Zhen’, LIU Shi-yin’, ZHANG Yong! , 
TANG Zhi-guang', WEI Jun-feng HUANG Dan-ni', ZHANG Sha-sha’ 
(1 Hunan University of Science and Technology , Xiangtan 411201 ,Hunan, China; 
2 School of Geomatics, Anhui University of Science and Technology , Huainan 232001 , Anhui, China; 


3 Institute of International Rivers and Eco-Security , Yunnan University , Kunming 650500 , Yunnan , China ) 


Abstract: The Karakoram glaciers are attracting more and more attentions for ‘ Karakoram Anomaly’ due to glac- 
ier surge or positive mass balance. To understand the mechanism of the glacier surge or positive mass balance, a 
huge debris-covered glacier — Yengisogat glacier was selected to be investigated in detail. Yengisogat glacier, the 
glacier of the largest area in China , developed on the northern slope of K2 located in the Karakoram Mountains, is a 
dendritic glacier covered area of 493 km? with four main tributaries and more than ten small glacier flows. TanDEM- 
X/TerraSAR-X data ( acquired in Feb. 2014) and SRTM-X dem ( acquired in Feb. 2000) were employed to calcu- 
late the glacier elevation change using differential interferometry method. Similar penetration depth was assumed for 
the same acquired season. The glacier surface elevation change and surface velocities from published articles were 
integrated to discuss the mass change characteristics. The results indicated that the glacier average surface elevation 
was decreased by 1. 68 +0.94 m between 2000 and 2014 , which means that the glacier has lost mass by - 0.12 + 
0.07 m * a '. The glacier surface elevation changes showed no uniformity. The southern main tributary was de- 
creased drastically. This was caused by rapid surface flow velocities which transfer the glacier mass from high alti- 
tude to low altitude ablation area. The southern main tributary occupied the original glacier trunk impeded the glac- 
ier mass transfer to low elevations which made the section elevation increased. The increased surface elevation in the 
tongue of glacier main trunk may be explained by glacier surge. Glacier surface elevation started to decrease along 
with the elevation increasing and then stayed stable , yet showing volatility at the same time. Surface elevation of thin 
debris-covered glacier in low altitude zones was prone to decrease than that of the exposed glacier ice. Glacier sur- 
face elevation was decreased more dramatically with the decreasing of glacier surface slope in the area where slope 
are less than 30 degrees. Surface elevation of glaciers with southern orientation was increased slightly while that with 
south-west orientation was decreased slightly, and the others were decreased drastically. The debris-covered 
Yengisogat Glacier experienced an overall surface elevation decline which indicated it had the mass loss although it 
had some local surface elevation increase due to the glacier surge. 


Key words:  Yengisogat Glacier; TDX/TSX ; D-InSAR ; surface elevation change; glacier surge 


